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("S , Two random- walk related problems which have been studied independently in the past, the ex- 

pected maximum of a random walker in one dimension and the flux to a spherical trap of particles 
undergoing discrete jumps in three dimensions, are shown to be closely related to each other and are 
studied using a unified approach as a solution to a Wiener-Hopf problem. For the flux problem, this 
work shows that a constant c = 0.297952f 9 which appeared in the context of the boundary extrap- 
olation length, and was previously found only numerically, can be derived analytically. The same 
^ , constant enters in higher-order corrections to the expected-maximum asymptotics. As a byproduct, 

. we also prove a new universal result in the context of the flux problem which is an analogue of the 

Sparre Andersen theorem proved in the context of the random walker's maximum. 
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, I. INTRODUCTION 

a . 

Random walks arise in an astounding variety of problems in physics as well as mathematics, computer science, 
etc., and great progress has been made in solving many of their deep and subtle properties. Two seemingly unrelated 
problems, the expected maximum of a random walker in one dimension undergoing jumps drawn from a uniform 
distribution, and the flux to a trap of particles undergoing random walks in three dimensions, have been studied 
independently over the last several years, and their solutions seem to involve a similar numerical constant 0.29795219, 
a coincidence that has not been noticed before. The similarity of these constants suggests that these two problems 
\Q I may be intimately related. For the first problem, the constant was first computed numerically by evaluating a rather 
Q\ • complicated double series expansion Q, Q and very recently, an exact closed form expression of the constant has been 
found 3] that is valid not just for the uniform jump distribution, but for any arbitrary continuous and symmetric 
jump distribution. For the second problem of flux to a spherical trap, the corresponding constant was computed only 
numerically Q . Therefore, finding the relation between the two problems raises the possibility that the flux problem 
can be solved analytically for the first time. Indeed, this is what we accomplish in this paper. 
The two problems we consider are: 
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Problem I: The asymptotic behavior of the expected maximum position of a discrete time random 
walker moving on a continuous line. The position x„ of the walker after n steps evolves for n > 1 via, 

O ■ 

O ' X n = X n -l + £ n (1) 

starting at xq = 0, where the step lengths £„'s are independent and identically distributed (i.i.d) random variables 
with zero mean and each drawn from the same probability distribution, Prob(£„ < x) — J_ f(y) dy, f(x) being 
a continuous and symmetric probability density normalized to unity. Let M n denote the positive maximum of the 
random walk up to n steps (see Fig. 1), 

M n =max(0,zi,x 2 ,-..,x n ). (2) 

We are interested in the asymptotic large- n behavior of the expected maximum E{M n ). This question arose some 
years ago in the context of a packing problem in two dimensions where n rectangles of variable sizes are packed 
in a semi-infinite strip of width one [3, 0- ^ was shown in Ref. that for the special case of the uniform jump 
distribution, f(x) = 1/2 for — 1 < x < 1 and f(x) = outside, for large n, 

E[M n ) = J — - 0.29795219028 • • ■ + 0(n" 1/2 ). (3) 
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FIG. 1: A typical configuration of a random walker in one dimension up to n steps, starting at at n = with M n denoting 
the maximum up to n steps. 



The leading s/n behavior is easy to understand and can be derived from the corresponding behavior of a continuous- 
time Brownian motion after a suitable rescaling 0. However, the leading finite-size correction term turns out to 
be a nontrivial constant — c with c = 0.29795219028 . . . that was computed in Ref. [3 by enumerating a somewhat 
intricate double series obtained after a lengthy calculation. Recently, it was shown |3| that for arbitrary continuous 
and symmetric jump distribution f(x) with a finite second moment a 2 — J Q x 2 f(x)dx, the expected maximum has 
a similar asymptotic behavior as in the uniform case, namely, 



E[M n ] = a 



c+Oin- 1 ' 2 ). 



(4) 



Moreover, an exact expression for the constant c was found 



7T 



dk , 



1 ~ /(fc) 

(T 2 fc 2 /2 



(5) 



where f(k) = J_ f(x) e lkx dx is the Fourier transform of f(x). In particular, for the uniform distribution, f(x) = 1/2 
for —1 < x < 1 and f(x) — outside, one has f(k) — sm(k)/k and one gets from Eq. © an exact expression, 



dk 

T 2 



In 



G 

k 2 



1 



sin k 



0.29795219028. . 



(6) 



Problem II: The calculation of flux to a spherical trap in three dimensions. Consider first the classic 
Smoluchowski problem |5j where point particles are initially distributed uniformly with density po outside a sphere of 
radius R in three dimensions. Each particle subsequently performs continuous-time Brownian motion with a diffusion 
constant D, independent of each other. One is interested in computing the flux of particles &(t) to the sphere at time 
t. This can be done by solving the diffusion equation outside the sphere with an absorbing boundary condition on 
the surface of the sphere and the result is well known d, IE 13 • One gets 



$(t) = 4irRDp 



1 



R 



y/irDt 



(7) 



valid for all t > 0. Also, as t — > oo, the density profile outside the sphere becomes time independent and has a simple 
form 



P( r ) = Pa [ 1 - — 
r 



(8) 
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FIG. 2: Independent Rayleigh flights in three dimensions in presence of a spherical trap with The discreteness of the jumps 
shows up effectively in a renormalized sphere with a smaller where c' I is the Milne extrapolation length 



for all r > R. Far from the sphere the density remains unchanged from its initial value po and as one approaches the 
surface of the sphere, the density vanishes. 

An interesting issue, first studied in Ref. Q, is how do the steady-state profile in Eq. (JHJ and correspondingly the 
expression of flux in Eq. J7J get modified when each of the point particles, instead of performing continuous-time 
Brownian diffusion, undergoes discrete 'Rayleigh flights', i.e. a particle jumps, at every discrete time step r, a fixed 
step length I whose direction is chosen arbitrarily in the three-dimensional space (see Fig. 2). In Ref. it was shown 
that the expression for the flux, at late times t = nr and for < I < 2R, now gets replaced by 



$(*) = ^k(R-c'1)D'pq 



1+ RJJ = + 0(t _3/ 2) 



(9) 



where D' = 1 2 /6t and c' is a constant whose numerical value was obtained by an iterative numerical solution of the 
density profile, with the result 

c' w 0.29795219 (10) 

The density profile p n (f) after n steps also gets modified rather drastically, as found numerically in Ref. Q. In 
particular, the steady-state density profile Poo(r) in the discrete problem turns out to be quite different from its 
continuous-time counterpart. While very far away from the sphere the steady-state density profile behaves as 

( R-c'l\ , x 

Poo(r»R) = p 1 (11) 



where c' is as in Eq. 110|) , the density actually tends to a nonzero constant as one approaches the surface of the sphere 
from outside 

Poo{r^R) = 0.408245^ (12) 
K 

where the constant 0.408245 was evaluated numerically in Ref. This is in stark contrast to the continuous-time 
Brownian case where the density vanishes on the surface of the sphere. 

The distance c' I is the 'Milne extrapolation length' 0, 0, and represents the distance inside the surface where 
the far steady steady-state solution in Eq. (|ll|l extrapolates to = 0. This steady-state density distribution implies 
the leading-order term of the flux given in Eq. (0. Thus at late times, the continuum formula for the flux in Eq. (J2J) 
still remains valid for the discrete jump case, but with an effectively smaller radius R— d 'I of the trap as in Eq. 0. 
Thus the effect on the flux due to the discrete nature of the jumps, at least at late times, is simply to renormalize the 
radius of the trap to a smaller value. 

Comparing Eqs. © and (|10fl one finds, rather amazingly, that the two constants c and c', in these two a priori 
unrelated problems, are identical at least up to 8 decimal places! This raises an interesting question: are they equal? 
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In this paper, we indeed prove that c = d . In the process, we also provide exact solutions to many other features of 
the flux problem that were observed numerically in Ref. [j]. For example, we will calculate exactly the steady-state 
density profile and will prove that indeed it approaches to a constant on the surface of the sphere as in Eq. (|12|l and 
the constant 0.408245 is actually l/v6- Our method consists in showing that both of these problems can be cast 
into the same Wiener-Hopf type problem involving an integral equation over half-space, albeit with different initial 
conditions. We then obtain explicit solutions to this Wiencr-Hopf problem with these two different initial conditions. 
The general solution turns out to be a product of two parts, one that explicitly depends on the initial condition and 
the other part which is a 'Green's function' that is independent of the initial condition. The constant c, given by the 
exact expression in Eq. ijjjj, is part of this 'initial condition independent' Green's function and hence it appears in 
both problems. 



II. MAXIMUM OF A RANDOM WALKER AS A WIENER-HOPF PROBLEM 



We consider a discrete time random walker hopping on a continuous line. The position x n of the walker evolves via 
Eq. JJJ, starting at xq — 0. The maximum M n , defined in Eq. is a random variable. Let q n {z) = Prob(Af„ < z) 
denote the cumulative distribution of the maximum and q' n {z) = dq n /dz its probability density with z € [0, oo]. Then, 
a simple integration by parts gives 



E[M n }= zq' n {z)dz = - z—(l-q n (z))dz = (l-q n (z))dz, 
Jo Jo az Jo 



(13) 



where we have used the normalization condition, q n {oo) = 1. Using translational invariance and the fact that the 
jump distribution is symmetric, one can also interpret q n (z) as the probability that a random walker, starting initially 
at position z > 0, stays positive up to step n. Then it is easy to write down, using the Markov property of the 
evolution in Eq. JTJ, the following recurrence relation Q, valid for all z > 0, 



q n (z)= / q n -i(z')f(z - z')dz', starting with q {z) = 1. (14) 
Jo 

The generating function q(z,s) — X)^Lo9™( z ) s ™ then satisfies an inhomogeneous Wiener-Hopf integral equation [jj 

poo 

q(z,s) = s / q(z',s)f(z- z')dz' + q (z), (15) 



where the inhomogeneous term qo{z) — 1 arises from the initial condition. We need to thus solve this integral equation 
to obtain the full probability distribution q n (z) of the maximum. The mean value E[M n ] can then be computed from 
Eq. Q. 



III. FLUX TO A TRAP AS A WIENER-HOPF PROBLEM 



We consider a sphere of radius R in three dimensions. Outside the sphere point particles are initially distributed 
with uniform density po. Particles subsequently perform independent Rayleigh flights, i.e. at every time step r, each 
particle, independently of others, jumps a fixed distance I in a direction chosen randomly. The object of interest is 
the flux at late times t to the sphere. For simplicity, we assume r = 1 (so that the jumps occur at integer steps) and 
also 1 = 1. Since the particles are independent, one can alternately think of a single particle whose probability density 
p n (r) at r after n steps evolves via the Markov equation Q 

Pn(rO = / W(f\f')p n ^(r')df' (16) 

where W(f\ r') = 8(\f— r'\ — l)/4ir is the jump probability density at each step from r' to r and the integral in Eq. 
(|l(j|| extends over the full three-dimensional space outside the sphere of radius R. The recursion relation in Eq. i|ltj|) 
starts with the initial condition, po(r) — pa for r > R and po(r) — for r < R. Since the initial condition is spherically 
symmetric, it is clear that from Eq. (|16|) that this symmetry will be maintained at all n, implying p n (r) = p n {f). 
Thanks to this spherical symmetry, the three-dimensional problem thus becomes a one-dimensional problem where 
one considers only the radial direction after integrating out the angular coordinates. Defining P n (r) = Aitr 2 p n {r) as 
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the probability that the particle is in the shell [r, r + dr] after n steps, it follows that P n {r) evolves via the recurrence 
equation 

P n (r) = / w(r\ r')P n -x(r')dr' starting with P (r) = 4nr 2 p 0(r - R) (17) 



where 0(x) is the Heaviside theta function and w(r\ r') is the jump probability density from radius r' to r, which can 
be easily calculated by integrating the kernel W(f\ r 1 ) over the angular coordinates pllTlj 

r dv 

w(r\ r )dr = -—j for \r — 1| < r < r +1 

= otherwise. (18) 

To simplify, one introduces a new quantity F n (r) = P n {r) / Ait p$r . The recursion relation for F n (r), upon substituting 
the explicit form of w(r\r') from Eq. (|18fl in Eq. I|17[). then simplifies 



1 



r+l 



F n (r) = - / F n -i(r f )dr f 

1 Jmax(Ji,r-l) 



F n - X {r')f{r-r')dr' (19) 



R 



where f(y) corresponds to the uniform probability density over the interval y G [—1, 1], i.e. f(y) = 1/2 if —1 < y < 1 
and f(y) = otherwise. The recursion in Eq. Ijl9(l now starts with the initial condition, Fo(r) = r0(r — R). One can 
simplify Eq. I|19l) further by introducing a shift, i.e. defining z = r — R and writing F n {r) = F n (z + R) = Q n (z), Eq. 
(fTT?fl becomes, for all z > 0, 

/>oo 

Qn(z)= Qn-i(z')f(z- z')dz', starting with Q (z) = R + z (20) 

Defining the generating function, Q(z, s) = J^^Lo Qn{z)s n , one obtains an identical Wiener-Hopf integral equation 
as in Eq. (fL5l) . 

poo 

Q{z,s) = s / Q(z',s)f(z-z')dz' + Q {z), (21) 



Jo 

the only difference is in the inhomogeneous term Q (z) = R + z that is set by the initial condition. One then needs 
to solve this integral equation to obtain Q n (z), from which one can read off the density profile at step n 

Pnir) = T^IT = ~ F n(r) = -Qn(r R). (22) 
The flux to the sphere at time step n can then be computed from the following relation Q 

$(n) = / dr' drw(r\r')P n - 1 (r') (23) 
Jr. Jo 

= irpo [ +1 [R 2 -(r' -l^Fn-^dr' (24) 
Jr. 

= irp f [2R(l~z)-(l-z) 2 ]Q n _ 1 (z)dz (25) 
Jo 

In going from Eq. Il2.'il) to l|24|l we have used the explicit form of w(r\r') in Eq. (JTSJl. 

We end this section with one remark. For a continuous-time Brownian motion it is quite standard 0, Il2| that, 
using the transformation p(r, t) = p F(r,t)/r, the 3-d diffusion equation for the density field p(r,t) can be reduced 
to a 1-d diffusion equation; the same trick naturally works for the 3-d Shrodinger equation as well. Based on this 
fact, it is natural that a similar transformation p n (r) — poF n (r)/r would also work for the discrete-time problem. 
However, the fact, that the reduced 1-d problem satisfies exactly the same integral equation (albeit with a different 
initial condition) with the same uniform kernel as the 1-d maximum displacement problem, is hard to guess apriori 
without the explicit calculation as presented here. 
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IV. WIENER-HOPF PROBLEM 



We have seen from the previous sections that the two a priori different problems (I) maximum of a random walker 
hoping on a line and (II) flux to a spherical trap in three dimensions can be both recast as the same Wiener-Hopf 
integral equation problem, albeit with different inhomogeneous terms arising due to the difference in the initial condi- 
tions of the two problems. The general mathematical problem then is to solve the following half-space inhomogeneous 
integral equation for z > 

P oo 

ij}(z,s) = s ip(z',s)f(z-z')dz' + J(z). (26) 
Jo 

where the inhomogeneous source term J(z) is different for the two problems 

J(z) = 1 for Problem I (27) 

= R + z for Problem II. (28) 

Even though in both of these problems the kernel f(z — z') corresponds to the uniform jump density, i.e. 

f(x) = - for - 1 < x < 1 

— otherwise, (29) 

it is useful to study the integral equation (|26[1 with a general continuous and symmetric kernel f(z) = f(—z) that is 
normalized f_ ao f(z)dz = 1 and has a finite second moment a 2 — J_ z 2 f(z)dz. 

The explicit solution ip(z, s) of Eq. I|26|) with different source terms as in Eqs. (|27|l and (|28fl will then provide the 
solutions to the two problems. In Problem I, ip(z,s) — q(z,s) = X)n=o Qn(z)s n provides the generating function for 
the cumulative distribution of the maximum of the random walk up to n steps. In Problem II, ip{z 1 s) — Q(z, s) = 
^2^=o Qn (z)s n gives the generating function for the density profile p n (r) — poQ n (z = r — R)/r of the particles outside 
the sphere of radius R in three dimensions. It turns out, as will be shown later explicitly, that the difference in the 
source term in Eqs. (|27|l and 128fl actually leads to completely different types of solutions to the integral equation 
(|26|l . In Problem I, the solution q n {z) depends explicitly on n even at late times, i.e. for large n, and does not have an 
n-independent stationary solution. Rather it has a scaling solution involving both z and n. In contrast, the solution 
Q n {z) in Problem II approaches an n-independent stationary solution. 



A. Explicit solution for exponential kernel 

Before providing the general solution for arbitrary continuous and symmetric kernel f(z), it is instructive to derive 
the explicit solutions with the two different source terms for a special kernel f(z) — exp[— |z|]/2. This will clearly 
bring out how the different source terms lead to different behavior of the same integral equation. The exponential 
kernel f(z) = exp[— \z\]/2 is special since one can recast the integral equation (|26[) into a differential equation by using 
the identity f"(z) — f(z) — S(z), where f"(z) = d 2 f/dz 2 . Differentiating Eq. I|2t)|) twice with respect to z and using 
the above identity, one gets for all z > 

| = (l- S )#, S ) + j"(z)-J(z). (30) 
Consider first Problem I where J(z) = 1. Then the most general solution of Eq. (|30|) is given by, 

M*, s ) = TT^ + Ms) e-^ z + Bx(s) e^ z , (31) 
(1 - s) 

where A\(s) and Bi(s) are two arbitrary z independent constants. Since the solution cannot diverge exponentially 
as z — > oo, one gets -Bi(s) = 0. We need to still determine the constant Ai(s). Here we use a method that is 
slightly different from that used in Ref. 0]. We substitute the solution in Eq. 1)31(1 into the integral equation l(26|l . 
Performing the integration explicitly, one finds that the solution in Eq. I|31l) satisfies the integral equation if and only 
if Ai(s) = —[1 — \/l ~ s \l (1 ~ s )- Thus, the full solution is given by Q 



m*, *) = E i^> n - tt-^t - i_ 9 s ?~^~ sz - ( 32 ) 

■n — n \ ) 
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One can then get the expected maximum from this explicit solution by an integration and an expansion in powers of 

, , 2 r(n + 3/2) Pfi , , 

E[M n ] = -= ■ ;/ -l~2W--l asn^co, 33 

V7r r(n + 1) V 1" 

which is of the same general form as in Eq. 10} with a = V2 and c = 1. In addition, it is also instructive to derive the 
solution q n (z) for large n by analysing its generating function in Eq. ll-il't in the vicinity of s = 1. Taking the limits 
s — > 1 and z — » oo but keeping z yl — s fixed, one can replace the the generating function by a Laplace transform, 
and inverting the Laplace transform one gets the scaling solution valid for large n 

q n (z) ~ erf (-±=) + -±= e^ 4 ", (34) 



where erf(z) = 2tt~ 1 / 2 J q 2 e~ u du is the error function. Note that the first term on the rhs of Eq. (|34|) corresponds 
to the continuum solution of the diffusion equation for a particle starting at z > and staying above an absorbing 
boundary at 0, which is also the same as the cumulative probability that the maximum of a continuous-time Brownian 
motion stays below z up to time t, provided one makes the standard correspondence a 2 n = 2Dt between the discrete 
step number n and the continuous time t, D being the diffusion constant for the Brownian motion. The second 
term on the rhs of Eq. I|34|) corresponds to the leading correction due to discrete jumps and indeed is responsible 
for the constant c in Eq. This can be seen by substituting the scaling solution in Eq. H34fl in the exact relation, 
E[M n ] = / °°(1 — q n {z))dz. Upon integrating, one recovers the large-n asymptotic solution in Eq. iffifl) and one sees 
explicitly that indeed the constant c = 1 in Eq. (|33|l arises from the integration of the second term in Eq. (f 341) . 

We now turn to Problem II where J(z) = R + z. Proceeding exactly as in the first case one finds that the explicit 
solution of the differential equation (|30|l is given by 

i; u ( z , s ) = ^±±+A 2 ( S )e-^ z , (35) 
1 — s 

where we have used the boundary condition that the solution cannot diverge exponentially as z — » oo. Substituting 
this solution in the integral equation l|26|> fixes the constant ^(s) and we get the full solution, 



M*> s ) = E ^{z)s n = ^ + (1 - R) 1 ,f\ S e-^ z . (36) 

The behavior of Q n (z) for large n can be derived by analysing the generating function near s — 1. In this case, one 
finds that for large n, 

Q n {z) ~ (z + 1) + OR - 1) erf ( -?=) + e~ z ^ 4n . (37) 



An 

Comparison with the asymptotic solution of Problem I in Eq. I|34() shows that in Problem II, the solution approaches 
an n-independent stationary solution as n — > oo 

Qoo(z) = z + 1 for all z > 0. (38) 

This solution is also independent of R; all terms containing R in Eq. (|37|l disappear in the long-time limit. The 
corrections to this stationary solution for large n have the scaling forms similar to Problem I. 

In fact, we will show in the next section that quite generically, i.e. for arbitrary continuous and symmetric kernel 
f(z), the solution Q n (z) for Problem II always approaches a stationary solution Qoo(z) which is, generically, a 
nontrivial function of z. However, for large z, we will show that this stationary solution has a rather simple asymptotic, 



Q 00 {z)~z + c' as z — » oo, (39) 

where the constant c' will be shown to be exactly equal to c in Eq. (J5J. In particular, for the uniform distribution 
f(z) given in Eq. (|29H where one can relate back to the original 3-d flux problem, we will show that indeed this same 
constant d — c appears as the extrapolation length in the expression for flux in Eq. 10 • In the particular example 
of the exponential kernel f(z) = exp[— |z|]/2, we see explicitly that indeed d = c = 1 by inspecting Eqs. (|38|l and 
(|33|l . This thus proves a special case of the general result d = c valid for arbitrary continuous and symmetric kernel 
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f(z). Note also that for this special case of the exponential kernel, the stationary solution Qoc(z) in Eq. I|38[l actually 
retains its asymptotic form in Eq. I|39|) all the way down to z = 0. This property, however, is rather special to the 
exponential kernel. For a generic continuous and symmetric kernel, Qoo{z) has a nontrivial form for small z as will 
be shown in the next section. 

Another quantity of interest, as we will see later in a more general context, is the transient behavior of Q n (0) for 
Problem II. It follows by subsituting z = in Eq. (|36|1 

00 1 D _ 1 

£Q n (0).s" = T - + ^ = =. (40) 

71=0 

Expanding the rhs of Eq. (|4U[I in powers of s one gets 

Q n (0) = l + (R-l)( 2n )-L^l + ^S. asn^oo. (41) 

\ Tl J Z a/ 7TTI 

We will see later that this transient behavior for the exponential kernel confirms, as a special case, the validity of a 
general result in Eq. I|9t)[l proved for arbitrary continuous and symmetric kernels. 



V. GENERAL SOLUTION TO THE WIENER-HOPF PROBLEM 



In this section, we present an explicit solution to the integral equation (|26[1 for the two different inhomogeneous 
terms in Eqs. (|27|l and (|28l) . Our result is valid for any arbitrary continuous and symmetric kernel f(z) that is 
normalized, /_ f{z)dz — 1 and with a finite second moment a 2 = J_ 00 z 2 j(z)dz. Our method relies on a beautiful 
general formalism developed by Ivanov 0] to deal with half-space problems in the context of photon scattering. Let 
us first summarize this formalism. Consider the integral equation l|26[l with an arbitrary source term J(z). There are 
three steps to obtain the solution. 

1. The first step is to define a Green's function G(z, zi, s) that satisfies the same integral equation but with a delta 
function source term, i.e. 

/•OC 

(42) 



(43) 



G(z, Zi,s) = s / G(z', zi f s)f{z — z')dz' + S(z — z\). 
Jo 

It is then easy to see that the solution of the inhomogeneous equation l|26() is given by 



ip(z,s) 



G(z, zi, s)J(z±)dzx. 



2. The next step is to determine the Green's function G(z, z\, s) that satisfies Eq. I|42fl. One first defines the double 
Laplace transform, 



G(A,A ljS ) 



dz e~ Az 



dzi e~ AlZl G(z,Zi,s). 



(44) 



Ivanov showed that this double Laplace transform can be determined in closed form by solving Eq. (|42|l and is 
given by a simple form 



G(A,Ai,s) 



( s > A) <j)(s, Ai 



A + Ai 

where the function <p(s, A) is the following Laplace transform 

/•oo 

0(s,A) = / dze- Xz G(z,0,s) 
Jo 

3. The third step is to obtain an explicit expression 0] for the function <f>(s, A) 

\ r ln(l-sf(k) 



(s, A) = exp 



A 2 + k 2 



dk 



(45) 



(46) 



(47) 



where f(k) = f(x) e lkx dx is the Fourier transform of the kernel /(a 
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Substituting the explicit expression for tfi(s, A) from Eq. i|47|) into Eq. (|45|l . one has an explicit expression for 
the double Laplace transform G(A, X±,s). By inverting this double transform, one can obtain the Green's function 
G(z, Zi, s), at least in principle. Subsequently, by performing the integral in Eq. H43JI one obtains the required solution 
i]j{z, s). In practice, however, these last two steps are difficult to carry out explicitly in general. However, for the two 
special source terms in Eqs. l{2"?|) and J35J), we show below that one can make progress. 



A. General solution for Problem I 



Consider first Problem I where J(z) = 1. Then, Eq. (|43[l gives 



G(z, zi, s)dz\. 



Let us define the Laplace transform 



tp I (X,s)= / ipi(z,s)e~ Az dz. 
Jo 



Taking the Laplace transform with respect to z in Eq. (|48|l we get 

/•oc />oo 

4>i(X,s)= dze- Xz G(z,z 1 ,s)dz 1 = G(\,0,s). 
Jo Jo 

Eqs. and JUJ then give 



^(A,s) = -0(s,A) <f>(s,0) 



(48) 



(49) 



(50) 



(51) 



where </>(s,A) is given in Eq. I|47|l . Let us first evaluate <j>(s, 0). Note that one cannot naively put A = in the 
expression in Eq. (|47|l since the integral multiplying A inside the exponential in Eq. I|47|l diverges as A — > 0. Hence one 
needs to extract the value of cf>(s, 0) carefully. To achieve this, an alternate expression for ^>(s,A) that was obtained 
in Ref. turns out to be useful. It was shown in Ref. that <f>(s, A) in Eq. (|47|l can also be written as 



( S ,A) 



exp 



dk 



In 



1 - sf(k) 



tt J Q A 2 + fc 2 \l-s + sa 2 k 2 /2 



(52) 



This representation is useful to derive the properties of <p(s, A) near A = 0. On the other hand, the representation in 
Eq. I)47|) is useful to extract the asymptotic behavior of (f>(s, A) for large A. Taking A — > limit in Eq. I)52[) one gets 



0M) 



VT 



which, when substituted in Eq. (|51|l gives 

MX,s) = y^s n q n (z)e~ Xz dz 

„_n JO 



n=0 



AVT 



■ <f>(s,X) 



(53) 



(54) 



where 4>(s, A) is defined in Eq. I|47|l and has also an alternative expression as in Eq. Q52JI. This result in Eq. I| 54ll goe s 
by the name of the Pollaczek-Spitzer formula which was originally derived using completely different methods |14lll5| . 
This result was subsequently utilized in Ref. [3j to extract the constant c in Eq. Q appearing as a subleading term 
for large n in the expected maximum E[M n ] of a random walker. 



B. General Solution for Problem II 

We now turn to Problem II where J(z) = R + z. We get from Eq. (|43(l 



V>n(z,s) 



(R + zi)G(z, zi, s)dz\. 



(55) 
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Taking the Laplace transform, ipu(X, s) = J ifin(z, s) e Xz dz gives 



iki(\,a) = I dze~ Xl 



[R + Zl ) G(z, z u s) dz, = R G(X, 0, s) - dG ^ X ^ s ) \ 



Eqs. (g5j) and gZJ then give 



■011 (A, s) 



X 



( S ,A) 



where 0(s, A) is given in Eq. (|47|) or alternately in Eq. I|52|) . Using the representation in Eq. I)52[| one gets 



d<f>(s,\ i) 
d\i 



Ui=0 



1 



(l-s)V2 ttVT 



00 dfc , 
ttt In 



1 - s />) 



s Jo ^ 2 \ 1 -s + scr 2 /c 2 /2 
Substituting this result in Eq. l(5T|) and using V>(s, 0) = 1/s/l — s from Eq. l(I)3"jl gives 



(56) 



(57) 



(58) 



^i(A,s) = i 



1 



1 - sf(k) 



— In I 

(1 - s) V 2 ' ttVT - ^ 7o fc 2 n \ 1 - a + sn 2 k 2 /2 



(s,X). (59) 



C. Analysis of the Results 

Thus the Laplace transform of the solution to the integral equation l|26l) for both problems respectively with J(z) = 1 
(Problem I) and J(z) = R + z (Problem II) can be written as a product of two functions 



$(A,s) = W(s, A) <f>(s, A) 



(60) 



where (p(s, A) is given in Eq. I|47|l or in Eq. I|52|) and is independent of the source term. The function W(s, A), however, 
depends explicitly on the source term, i.e. on the initial conditions of the original recursion relations in Eqs. i|14fl and 
(I2U|1 and has different expressions for the two problems. While for Problem I it is rather simple 



Wi(a,X) 



XVT 



for Problem II it has a more complicated expression 



Wu(s,X) = - 



1 



pi 

VT~s V x) ' (1 - s) V 2 ' k^TS 7 k 2 



1 



00 dk , 

TTT hi 



1 - sf(k) 
1- s+ sa 2 k 2 /2 



(61) 



(62) 



Since the function <fi(s, X) appears in both problems, it is useful to first list its asymptotic properties for small and 
large A. For small A, it is convenient to use the representation in Eq. (|52|) . On the other hand, for large A, the 
representation in Eq. I|47|) turns out to be more convenient. We find 



(s,X) 



1 



(3(s) 



1 



A 



-a(s)X + 0(X 2 ) asA^O 
+ 0(A~ 2 ) asA-^oo 



where the two functions a(s) and (3(s) are given by 



a(s) 



a 



m = - 



(1-8) 



1 



1 - Sf(k) 



TV 



dk In 



o 



— 1 1 

nVT^Jo k 2 n \l- S + sa 2 k 2 /2 



1 - sf(k) 



(63) 
(64) 

(65) 
(66) 



Now we are ready to obtain the constants c in Problem I and c' in Problem II and show that indeed they are same, 
as they both emerge from the properties of the function 0(s, A) that is common to both the problems. We consider 
the two problems separately. 
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Expected maximum in Problem I: Consider first Problem I. Consider the Laplace transform of the distribution 
of the maximum, E[e~ XAIn ] = e~ Xz q' n (z) dz = A J °° e~ Xz q n (z) dz. It follows from Eq. lt5^|) 



-AM. 



n=0 



(s,X). 



(67) 



We expand both sides in A for small A and use Eq. (|63|) for the expansion of the rhs of Eq. (|67|1 . Comparing the term 
linear in A one gets 



J2 s n E[M n 



71=0 



(68) 



where a(s) is given in Eq. (|65|l . To extract the large n behavior of E[M n ], we need to analyse the generating function 
in Eq. (|68|) near s — > 1. Expanding the rhs of Eq. (|68|) near s — 1 one gets 



V s n E[M n ] ~ = °" 
^ V2(l- 

n— v 

Inverting one gets the large n result Q 





"(!-/»" 




CT 2 fc 2 /2 



+ 0((1- S )-V2). 



2?[M n ] =C7 A /— -c+0(ti- 1 / 2 ), 



(69) 



(70) 



where the constant c is given by the integral in Eq. 

Steady-state density profile in Problem II: Turning now to Problem II, we first show that for large n, Q n (z) 
approaches a stationary solution Qoo{z). From the definition, we have 

5>"/ Qn{z)e~ Xz dz = ^n(A,a), (71) 

n=Q Jo 

where -0n(A, s) is given in Eq. Ij59|) . Thus, if Q n (z) — ► Qoo(^) as n — > oo, then the lhs of Eq. Ij71(l will behave as 

Vs" g„(2)e- Az ^~ / Q 00 (z)e- Xz dz as s^l. (72) 

We now expand the rhs of Eq. I|71|l near s = 1 using the explicit expression of tpn(X, s) in Eq. I|59fl . We find from Eq. 
(|59|l that as s — > 1, the leading order term behaves as 



■4>n(A,s) 



a(f>(l,X) 



(73) 



v^A(l - s) 

Comparing Eqs. I|72l) and (|73[) gives the exact Laplace transform of the stationary solution 

'QMe- x *dz=^l, 
A v 2 

where 0(1, A) can be obtained by putting s — 1 either in Eq. (|47|) or alternately in Eq. I|52|) . Both expressions are 
equivalent and one gets 



(74) 



3 (z)e~ Az dz 



Aa/2 
1 
A 



exp 



A 



dk 



^exp 



ir J A 2 + fc 2 
A f°° dk 



In (l - f(k) 
A 2 + fc 2 ~~ I cr 2 fc 2 /2 i 



In 



(75) 
(76) 



As a check, one can verify that for the exponential kernel f(x) = exp[— \x |]/2 so that f(k) = 1/(1 + fc 2 ) and a = y/2, 
Eq. gives 



(z)e- Xz dz= ^ + i 



(77) 
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which, when inverted, gives Qoc(z) = z + 1, in agreement with Eq. (|38|) . 

For a general kernel it is difficult to invert the Laplace transform in Eqs. I|75|l or H76|) . However, one can easily 
extract the asymptotic behavior for large and small z. Consider first the large z behavior. In this case, the expression 
in Eq. (|76[) is more convenient. Expanding the rhs of Eq. I|76[) for small A we get 



1 



QMe-^dzc--— I -In 



Inverting the Laplace transform, one gets 



where 



Qoo(z) ~ z + c' 



dk 



k 2 



as z — > oo 



<J 2 k 2 /2 



dk 



In 



1 - />) 
a 2 k 2 /2 



(78) 



(79) 



(80) 



thus proving one of the main results of this paper. In particular, for the uniform kernel in Eq. i|29[l such that 
f(k) = sin(fc)/fc, we get 



c = 0.29795219028.. 



(81) 



The asymptotic exact result Qoo(z) — > z + c' as z — > oo with c' = 0.29795219028.., when substituted in Eq. l|22|l in the 
steady-state n — > oo limit, provides the exact steady-state density profile in the original 3-d flux problem at distance 
r » R, 



Poo (r » R) = ^Qoo(r - R) = p ( 1 - 

r \ r 



(82) 



in perfect agreement with the numerically observed Q| behavior in Eq. i|ll|) (note that we have set I = 1 in Eq. I|ll|)l. 

Similarly, one can work out the small z behavior of Qoo{z) by analysing the large A behavior of the Laplace 
transform. In this case, the expression in Eq. {75} is more convenient. Expanding Eq. (|75|l for large A one finds 



Qoo{z)e- Az dz = -j= 



dk\n[l-f(k))+0(\- 3 ) 



Inverting the Laplace transform, we get 



where & is a new constant given by 



Qoo(z) ~ — + bz 



as z — ► 0, 



dk ln(l - f(k) 



For the exponential case, f(k) = l/(k 2 + 1) and a — y/2 we get from Eq. I|85|) 



/ dk In 


'(fc 2 + l)" 


= 1 


/o 


fc 2 





(83) 



(84) 



(85) 



(86) 



Thus from Eq. I|84l) . for the exponential kernel, we get Qoo(z) = z + 1 as z — ► 0, in agreement with Eq. I|38|) . For the 
uniform case, f(k) = sin(k)/k and a = l/y/3, we get 



b= -- 





^ sin(fc) 


/ dk In 


/o 


k 



= 0.653857. 



(87) 



Thus, for the uniform kernel, the small z behavior of Qoo(z) from Eq. (|84|) reads as 



Q oc (z)~ ^ + (0.653857... )z as z -> 
V6 
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This is in per fect agreement with the numerical results obtained from the work of Ref. [4|, Qoo(z) ~ 0. 408245 + 0. 6538z 
as z — ► [Ig- Substituting the result from Eq. (|88|l in Eq. (|22[) in the limit n — > oo, we find the steady-state density 
profile in the original 3-d flux problem near the surface r —> R 

Poo(r -^R) = ^Qoo{r - R) w ^ + 0.653857 ■ • • (r - . (89) 

In particular, on the surface, the steady-state density approaches to a constant value 

Poo (R) = .Le°, (90) 

thus proving the numerically observed Q relation in Eq. \12\ and identifying the constant 0.408245 as l/VQ (note 
that we have set I = 1 in Eq. iJHJ). 

Transient behavior in Problem II: Another quantity that was investigated numerically in Ref. j4| is the transient 
behavior of Q n (0) for large n, where it was found that for the uniform kernel in Eq. I|29l) 

Q„(0) fa 0.408245 + [-= - 0.17 ) n~ 1/2 (91) 



Our exact solution reproduces this result also. To see this, we investigate our explicit solution in Eq. (|59|l for a general 
continuous and symmetric kernel. By taking A — > oo limit in Eq. 171|) we see that the lhs behaves as 

J2 sU / Qn{z)e- Xz dz ~ -Vg„(0)s" asA^oo. (92) 

n=0 "'O A n=0 

On the other hand, expanding the exact expression of s) in Eq. (|59|l for large A, we find that the rhs of Eq. I|71|) 

behaves as 



V>n(A,s) ~ - 



, , R 

a{s) + 7T= 



as A — > oo (93) 



where a(s) is given in Eq. (HI). Equating Eqs. and §53J| and using the expression for a(s) from Eq. (|65|l we get 
an exact expression for the generating function 

„ ° [s 1 f°° dk ( l-sf(k) \ R 

g<3„(0K^_ v /- + - 7 =y () F h ^_ t + ;u 2/2 j + _ . ,94) 

Now, analysing the behavior near s = 1 of the rhs of Eq. I|94|). one can get the leading asymptotic behavior of Q n (0) 
for large n. One gets the leading behavior near s = 1, 

where c is the same constant as in Eq. (|5j). Inverting, one obtains an exact asymptotic result for large n 

V2 V7T 

In particular, for the uniform kernel /(/c) = sin(fc)/fc with a = l/v3 and c = 0.29795219028 . . . from Eq. I|81|) we get 

^ 1 + (^- 0.29795219028...) ^ = Q + / R _ \ ^ 

V6 Vtt VV 71 " / 

in excellent agreement with the numerical results stated in Eq. I|91|) . 



14 



VI. SPARRE ANDERSEN THEOREM FOR PROBLEM I AND ITS ANALOGUE FOR PROBLEM II 

The recursion relation in Eq. (|14fl satisfied by the q n (z)'s in Problem I has an explicit solution given in Eq. 154fl . 
From this explicit solution, one can easily extract q n (fi), the probability that a random walker starting at the origin 
stays positive (or negative) up to n steps. Indeed, the integral J °° q n (z)e~ Xz dz — > g n (0)/A in the A — > oo limit. On 
the other hand, the rhs of Eq. (|51Jl tends to l/A\/l — s as A — > oo since cj>(s, oo) = 1. Matching the lhs and the rhs 

gives, 

oo - 

]>>„(0) S " = -=. (98) 
Expanding the rhs of Eq. (|98(l in powers of s and identifying the coefficient of s n on both sides, one gets for all n 

^(0)=( 2 ;)^. (99) 

The amazing fact is that the result in Eq. (|99|l is universal for all n, i.e. it does not depend on the density function 
f(z) as long as f(z) is continuous and symmetric. This, in fact, is the celebrated Sparre Andersen theorem which was 
originally proved using combinatorial methods |17| and has since been reproduced by various other methods [TR Il8j . 
In particular, one notes from Eq. I|99|l that in the limit of large n, 

lim \fn q n (0) = —j= = a universal constant (100) 

n— too yjn 

A natural question is if there is an analogue of the universality a la Sparre Andersen theorem for Problem II. The 
recursion relation for Problem II in Eq. (|2U[) is identical to that of Problem I in Eq. (|14f> . except that the initial 
condition Qq(z) = R+ z is different from that in Eq. (|14fl . The question is whether Q„(0) still remains universal with 
this different initial condition. Indeed, the answer to this question is evident from our exact result in Eq. H94J) . It is 
clear from Eq. I|94(l that unlike in Problem I, Q n (0) in Problem II is not universal for all n as it depends explicitly on 
the density function f(k). However, one recovers universality (up to a constant scale factor a ) asymptotically, i.e. in 
the limit of large n. Indeed, it follows from Eq. H96|l 

lim = —= = a universal constant (101) 

n—>oo a y/2 

Thus Q„(0), scaled by a, approaches a universal constant 1/V2 as n — > oo, independently of the density function 
f(z) as long as f{z) is continuous, symmetric and has a finite second moment a 1 = j_ z 2 f(z) dz. The result in 
Eq. i|l(Jl[l for Problem II can be thought of as an asymptotic analogue of the Sparre Andersen result in Eq. (|1(J0[I for 
Problem I. 



VII. CONCLUSIONS 



We have shown that the two apparently different discrete-time random walk problems, the expected maximum in 
one dimension and the three-dimensional flux to a trap, are intimately related in that they end up satisfying the 
same recursion relation, but with different initial conditions. We confirmed that the constant c' in the flux problem is 
identical to the constant c in the expected-maximum problem, and thus provide an exact derivation for this constant 
which was previously found only numerically, fourteen years ago. We also find the surprisingly simple result that 
the steady-state density in the flux problem reaches a constant at the boundary of the sphere, p 00 (R) = pol/(\/6R) 
or equivalently Qoo(0) — 1/V6, and find explicitly both the asymptotic time-dependent approach to that value, Eq. 
(19 711 . and the slope of the density at the wall at z = in the steady state, Eq. (|88|l . 

For the flux problem, c' represents the extrapolation length inside the boundary where the steady-state solution 
far from the sphere Poo{r) = pa{r — R + c'£)/r goes to zero. That is, the solution far from the sphere assumes the 
form of the solution to the diffusion equation, but with the effective boundary somewhat inside the actual boundary. 
Putting this expression into the formula for the flux, $ = 4Trr 2 D(d/dr)p oc (r), yields the leading term in Eq. @ (with 
1 = 1). We evaluate $ for large r where the solution is valid, and the result is naturally independent of r because it 
is in steady state. 

In this paper, we have presented explicitly the steady-state density profile for the 3-d discrete flux problem. As 
mentioned above, this solution suffices to predict immediately the leading behavior of the flux $(£) as t — > oo in Eq. 
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To obtain the subleading time-dependent term given in Eq. J^J, it is necessary to study the large time asymptotic 
behavior of the integral in Eq. i|25[l. Indeed, using the asymptotic solution for Q n (z) derived here, it is possible 
to calculate this subleading term explicitly; we have not presented this calculation here and it will be published 
elsewhere [l9j| . It is worth remarking that in the result found for the flux in Eq. J^J, the combination R — c'l enters in 
the time-dependent correction in the same form as in the steady-state term — that is, the extrapolation-length idea 
also figures in the asymptotically large-time behavior of the system. 

We have considered the three-dimensional flux problem here, in which particles undergo a Rayleigh flight of unit step 
length. The resulting equation, Eq. I|2U|) . can also be interpreted as representing a one-dimensional flux problem where 
the particles undergo a uniformly distributed jump, with an adsorbing boundary at z — 0. For the one-dimensional 
interpretation, however, the initial condition is uniform rather than linear as in l|28[) . (Call this Problem II'). Then 
Problems I and IF become identical, with E[M n ] of the former, Eq. (|13|) . corresponding to the total integrated flux 
up to time n of the latter. For a uniform jump distribution, Eq. @ thus gives the accumulated flux up to time n. In 
Ref. Q, the quantity a^ 2 \n) of Eq. (10a) is equal to six times the flux at time n for precisely this one-dimensional 
problem. Summing the data for a^ 2 \n)/d> up to time n (some of which data is presented in Table I of that paper), 
we indeed find that Eq. © (including c) is satisfied, with the next correction approximately equal to 0.0921n -1 / 2 . In 
fact, the latter agrees with the prediction (l/5)-\/2/37m = 0.09213177. . .rT 1 ! 2 given in Ref. @ and this prediction 
also agrees with the conjectured behavior of (n) given in Eq. (21c) of Ref. 0. Thus, we have verified an additional 
conjecture of Ref. 0. 

On the mathematical side, we have discussed explicit solutions to the Wiener-Hopf type integral equations for 
continuous and symmetric kernel with a finite second moment a 2 . This is because both the physics problems that we 
were interested in this paper have kernels that satisfy these properties. Mathematically it is interesting to ask how 
the solutions would change if the kernel is asymmetric or for example, does not have a finite second moment. This 
later problem with diverging second moment corresponds to discrete-time Levy flights and at least for the expected 
maximum problem (Problem I), exact results using similar techniques have recently been obtained On the other 
hand, while general solutions to the Wiener-Hopf type integral equations with asymmetric kernels can, in principle, 
be obtained [l3j, they are mostly not explicit. Finding explicit solutions with asymmetric kernels remains a hard and 
challenging problem. 

Finally, our solution of the recursion relation in Eq. Ijl4|l with a constant initial condition includes, as a special 
case, a simple derivation of the Sparre Andersen theorem that states that q n (0), he. the probability that starting at 
a random walker's path stays above (or below) up to n steps, is independent of the jump density function as long 
as it is continuous and symmetric. However, we have proved that for the same recursion relation, but with a linear 
initial condition as in Eq. (|20|l . the universality in Q rl (0) holds only asymptotically for large n up to a constant factor, 
i.e. the ratio Qn(0)/cr (where a is the standard deviation associated with the jump density function) approaches to 
a universal constant 1/V2, irrespective of the details of the jump density function as long as it is continuous and 
symmetric. This paper provides a rigorous proof of this new theorem. 
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